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Context. The seismic data obtained by CoRoT for the star HD 49933 enable us for the first time to measure directly 
the amplitudes and linewidths of solar-like oscillations for a star other than the Sun. From those measurements it is 
(— I ' possible, as was done for the Sun, to constrain models of the excitation of acoustic modes by turbulent convection. 

Aims. We compare a stochastic excitation model described in Paper I with the asteroseismology data for HD 49933, a 
I , star that is rather metal poor and significantly hotter than the Sun. 

O ' Methods. Using the seismic determinations of the mode linewidths detected by CoRoT for HD 49933 and the theoretical 

^ ' mode excitation rates computed in Paper I for the specific case of HD 49933, we derive the expected surface velocity 

t/3 , amplitudes of the acoustic modes detected in HD 49933. Using a calibrated quasi-adiabatic approximation relating the 

. mode amplitudes in intensity to those in velocity, we derive the expected values of the mode amplitude in intensity. 

Results. Except at rather high frequency, our amplitude calculations are within 1-a error bars of the mode surface 
velocity spectrum derived with the HARPS spectrograph. The same is found with respect to the mode amplitudes 
in intensity derived for HD 49933 from the CoRoT data. On the other hand, at high frequency {v > 1.9 mHz), our 
calculations depart significantly from the CoRoT and HARPS measurements. We show that assuming a solar metal 
abundance rather than the actual metal abundance of the star would result in a larger discrepancy with the seismic 
data. Furthermore, we present calculations which assume the "new" solar chemical mixture to be in better agreement 
with the seismic data than those that assumed the "old" solar chemical mixture. 

Conclusions. These results validate in the case of a star significantly hotter than the Sun and a Cen A the main 
■ assumptions in the model of stochastic excitation. However, the discrepancies seen at high frequency highlight some 

deficiencies of the modelling, whose origin remains to be understood. We also show that it is important to take the 
^\ [ surface metal abundance of the solar-like pulsators into account. 
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1. Introduction ing into account the measured iron-to-hydrogen abundance 

Cd ' ([Fe/H]) of HD 49993 ([Fe/H]= -0.37), we have derived 

• The amplitudes of solar-Hke oscillations result from a bal- the theoretical values of the mode excitation rates ? ex- 

ance between excitation and damping. The mode linewidths pcctcd for this star. The resulting value of V is found to be 

are directly related to the mode damping rates. Once we can ^bout two times smaller than for a model with the same 

measure the mode linewidths, we can derive the theoretical gravity and effective temperature, but with a solar metal 

value of the mode amplitudes from theoretical calculations abundance (i e [Fe/H]= 0) 

of the mode excitation rates, which in turn can be compared xhc s tar HD 49933 was fi rst observed in Doppler ve- 
to the available seismic constraints. This comparison allows lo^ity by iMosser et al.] (l2505li with the HARPS spectro- 
us to test the model of stochastic mode excitation investi- g^aph. More recently, this star has been observed twice 
gated in a companion paper (|Samadi et al.||2009|, hereafter CoRoT. A first time this was done continuously during 
Paper I). about 61 days (initial run, IR) and a second time contin- 
As shown in Paper I, a moderate deficit of the sur- uously during about 137 days (first long run in the ccn- 
face metal abundance results in a significant decrease of ter directio n, LRcOl). The comb ined seismic analysis of 
the mode driving by turbulent convection. Indeed, by tak- these data (jBenomar et al.l |2009[) has provided the mode 
linewidths as well as the amplitudes of the modes in inten- 



* The CoRoT space mission, launched on December 27 2006, sity. Then, using mode linewidths obtained for HD 49933 

has been developped and is operated by ONES, with the contri- with the CoRoT data and the theoretical mode excitation 

bution of Austria, Belgium, Brasil, ESA, Germany and Spain. rates (obtained in Paper I), we derive the expected values 

Correspondence to: Reza.Samadi@obspm.fr of the mode surface velocity amplitudes. We next compare 
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these v alues with the mode v elocity spectrum derived fol- 
lowing iKieldsen et al.l (l2005f) with seismic data from the 
HARPS spectrograph (jMosser et al . 2005). 

Mode amplitudes in terms of luminosity fluctuations 
have also been derived from the CoRoT data for 17 radial 
orders. These data provide us with not only a constraint 
on the maximum of the mode amplitude but also with the 
frequency dependence. The relative luminosity amplitudes 
SL/L are linearly related to the velocity amplitudes. This 
ratio is determined by the solution of the non-adiabatic pul- 
sation equation s and is independent of the stochastic excita- 
tion model fsee lHoudek et aLlllQQOt ). Such a non-adiabatic 
calculation requires us to take into account, not only the 
radiative damping, but also the coupling between the pulsa- 
tion and the turbulent convection. However, there are cur- 
rently very significant uncertainties concer ning the mode l- 
ing of this coupling (for a recent review see lHoudekll2008l ). 
We relate further for the sake of simplicity the mode lumi- 
nosity amplitudes to computed mode velocity amplitudes 
by ass uming adiabatic oscillations as IKieldsen fc Bedding] 
(jl995f ). Such a relation is calibrated in order to reproduce 
the helioseismic data. 

The comparison between theoretical values of the mode 
amplitudes (both in terms of surface velocity and intensity) 
constitutes a test of the stochastic excitation model with 
a star significantly different from the Sun and a Cen A. 
In addition it is also possible to test the validity of the 
calibrated quasi-adiabatic relation, since both mode ampli- 
tudes, in terms of surface velocity and intensity, are avail- 
able for this star. 

This paper is organized as follows: We describe in Sect. [2] 
the way mode amplitudes in terms of surface velocity Vs 
are derived from the theoretical values of T and from the 
measured mode linewidths (F) . Then, we compare the the- 
oretical values of the mode surface velocity with the seismic 
constraint obtained from HARPS observations. We describe 
in Sect. [3] the way mode amplitudes in terms of intensity 
fluctuations SL/L are derived from theoretical values of 
and compare 6L/L with the seismic constraints obtained 
from the CoRoT observations. Finally, Sects. [4] and [5] are 
dedicated to a discussion and conclusion respectively. 

2. Surface velocity mode amplitude 

2.1. Derivation of the surface velocity mode amplitude 

The intrinsic rms mode surface velocity Vg is related to 
the mode exitation rat e J'(t^) and the mod e linewidth F(;^) 
according to (see, e.g.. lBaudin et al.|l2005f ): 



27r F 



(1) 



where CP is the mode excitation rate derived as described in 
Paper I, F is the mode full width at half maximum (in u)^ 
V = LOoscl'2'T^ the mode frequency and is the mode mass 
defined as: 

M,. = (2) 

where / is the mode inertia (see Eq. (2) of Paper I), ^j- the 
radial mode eigendisplacement, = R + h the layer in the 
atmosphere where the mode is measured in radial velocity. 
R the radius at the photosphere (i.e. at T = Toff) and h the 
height above the photosphere. 



In Sect. [272] we will compare e stimated values of Vs with 
the seismic constraint obtained bv lMosser et all (120051 ) with 
the HARPS spectrograph. We therefore need to estimate 
Vs at the layer h where the HARPS spectrograph is the 
most sensitive to the mode displacement. As discussed by 
I Samadi et al.l (|2008aD . the seismic measurements obtained 
with HARPS spectrograph are likely to arise from the opti- 
cal depth r 500 nm — 0.013, which corresponds to the depth 
where the potassium (K) spectral line is formed. We then 
compute the mode mass at the layer h associated with the 
optic al depth r 500 nm (|Christensen-Dalsgaard fc Goughl 
119821). For the model with [Fe/H]= (resp. [Fe/H]= -1) 
this optical depth corresponds to h ~ 390 km (resp. h ~ 
350 km). 

For the mode linewidth F we use the seismic measure- 
ment obtained from the seismi c analy sis of the CoRoT data 
performed by iBenomar et al.l l|2009f l. This seismic analy- 
sis combined the two CoRoT runs available for HD 49933. 
Two different approaches were considered in this analysis: 
one based on the maximum likelihood estimator and the 
second one using the Bayesian approach coupled with a 
Markov Chains Monte Carlo algorithm. The Bayesian ap- 
proach remains in general more reliable even in low signal- 
to-noise conditions. Nevertheless, in terms of mode ampli- 
tudes, mode heights and mode linewidths, both methods 
agree within 1-a. We will consider here the seismic param- 
eters and associated error bars obtained on the basis of the 
Bayesian approach. 

2.2. Comparison with the HARPS measurements 

The seismic analys is in velocity has been performed by 
iMosser et"all (|2005l ) using data from the HARPS spectro- 
graph. The quality of these data is too poor to perform 
a direct comparison between the observed spectrum and 
the calculated amplitude spectrum {vs, Eq. ([T])). Indeed, 
the observed spectrum is highly affected by the day aliases. 
Furthermore, the quality of the data does not allow to iso- 
late individual modes, in particular modes of a different an- 
gular degree {£). A consequence is that energies of modes 
which are close in frequency are mixed. 

In order to measure the oscillation amplitude in a 
way that is independent of these effects, we have followed 
the method in t roduc ed by IKieldsen et al.. (20051 . see also 
iKjeldsen et al.l (|20G8f)). This method consists in deriving 
the oscillation amplitudes from the oscillation power den- 
sity spectrum smoothed over typically four times the large 
separation (i.e. four radial orders). Next, we multiply this 
smoothed spectrum by a coefficient in order to convert the 
apparent amplitudes into intrinsic amplitudes. This coef- 
ficient takes into account the spatial respo nse function of 
the angular degrees i =0, 1, 2 and 3 (see IKieldsen et all 
200§). We have checked that the sensitivity of the visi- 
bility factor with the limb-darkening law is significantly 
smaller in comparis on with the error associated with the 
iMosser et al.l (j2005l ) seismic measurements. The amplitude 
spectrum wharps derived foUowing Kieldsen et al. (2005) is 
shown in Fig. |TJ The 1-a error bar associated with each val- 
ues of wharps is constant and equal to AynARPS = 7 cm/s. 

The maximu m of t;HAR.PS r e aches Vmax = 50.2 ± 7 cm/s. 
By comparison, IMosser et al.l (j2005l) found a maximum of 
40 ± 10 cm/s, which once converted into intrinsic ampli- 
tude represents a maximum of 42 ± 10 cm/s. The differ- 
ence between the two values is within the 1-a error bars. 
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The difTcrcnt value found by iMosser et aLl ()2005[ ) can be 
explained by the way the maximum of the m ode ampli- 
tude was derived. Indeed, iMosser et al.l (|200S1 have con- 
structed synthetic time series based on a theoretical low de- 
gree p-modes eigenfrequency patte rn and theoretical mode 
lines widths (jHoudek et al.l [l999l ). The maximum ampli- 
tudes were assumed to follow a Gaussian distribution in 
frequency. Using a Monte-Carlo approach, the maximum 
amplitude was then determined in order to obtain com- 
parable energy per frequency bin in the synthetic and ob- 
served spectra. On the other ha nd, except for the mo de re- 
sponse function, the method bv lKieldsen et all (|2005[ ) does 
not impose a priori constraints concerning the modes. This 
method can t hen be considered to be more reliable than 
the method bv IMosser et all (|2005[ ). 

We compare in Fig.[T]wHARPS with the calculated mode 
surface velocity Vg (Eq. However, in order to have a 
consistent comparison, we have smoothed Vg quadratically 
over four radial orders. We note Avs the 1-a error bars as- 
sociated with Vs- They are derived from AT, the 1-cr error 
bars associated with T. As pointed out in Paper I, the un- 
certainty related to our knowledge of the metal abundance 
Z for HD 49933 results in an uncertainty about the deter- 
mination of However, in terms of amplitude, this uncer- 
tainty is of the order of 5 % ; this is negligible compared to 
the uncertainty that arises from AT (ranging between 25 % 
to 50 % in terms of amplitude). 

The difference between computed values and observa- 
tions is shown in the bottom panel of Fig.[T] This difference 
must be compared with a^, the 1-a interval resulting from 
the errors associated with Vg and this in turn associated 
with tiHARPS, that is ay = \/ Aw^ -|- Aw^^ppg. As seen in 
Fig. d except at high frequency {y^\.% mHz), the theo- 
retical Vs lie well in the l-cr„ domain. However, there is a 
clear disagreement at high frequencies where the computed 
mode surface velocities overestimate the observations. This 
disagreement is attributed to the assumptions in the theo- 
retical model of stochastic excitation (see Sect. 14. 5| ). 

Assuming the 3D model with the solar abundance re- 
sults in significantly larger Vs- In that case the differences 
between computed Vs and the seismic constraint are in gen- 
eral larger than 2-cr^. This shows that ignoring the metal 
abundance of HD 49933 would result in a larger discrepancy 
between Vs and wharps- 



3. Amplitudes of mode in intensity 

3.1. Derivation of mode amplitudes in intensity 

Fluctuations of the luminosity L due to variations of the 
stellar radius can be neglected since we are looking at high 
n order modes; accordingly the bolometric mode intensity 
fluctuations SL are mainly due to variations of the effective 
temperature, that is: 



L 



(3) 



As in lKieldsen fc Beddind (|1995[) . we now assume that STcs 
is proportional to the variation of the temperature induced 
by the modes at the photosphere (i.e. at T = Tcff). This 
assumption is discussed in Sect. 14.31 Assuming further low 
degree £ and adiabatic oscillations, one can derive a relation 
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Fig. 1. Top: Intrinsic mode surface velocity as a function of 
the mode frequency {v). The filled circles connected by the 
thick solid line correspond to the mode surface velocity (vs) 
derived for HD 49933 according to Eq. ([1]), where the mode 
excitation rates CP are derived as explained in Paper I and 
the mod e linewidths and their associated error bars are de- 
rived bv lBenomar et al.l (|2009D from the CoRoT data. The 
thick dashed line corresponds to the mode velocity asso- 
ciated with the model with [Fe/H]= 0. The thick and red 
solid line corresponds to the amplitude spectrum derived 
from the seismic observations obtained with the HARPS 
spectrograph (see text). The dotted line corresponds to the 
1-a domain associated with this measurement. Bottom: 
Differences betwee n Vs and iiharps- The 1-a error bars cor- 
respond to av = + Aw^^j^pg (see text). 



between STcs/Tcs and the radial mode velocity v that is: 



ST, 



eft 



off 



(Fa - 1) 



_d6 

dr 



(4) 



where = Vad Fi + 1, Vad is the adiabatic temperature 

' d In P„ 



gradient, Fi = 



din p 



the radial mode eigendisplace- 



ment, and v the mode velocity at the photosphere. Finally, 
according to Eqs. ^ and one has: 



1) 



1 dgr 



(5) 
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where v is computed using Eq. ([T]) with h~0 (the photo- 
sphere), that is: 



associated with 6L/L. Accordingly, we have A{SL/L) = 



2n MoT 



(6) 



where Mq is the mode mass evaluated at the photosphere 
(h^O). 

In Eq. ([5]), /3 is a free parameter introduced so that 
Eq. ^ gives, in the case of the solar p modes, the cor- 
rect maximum in SL/L. Indeed, Eq. ([5|) applied to the 
case of the solar p modes, overestimates by ~ 10 times the 
mode amplitudes in intensity. This important discrepancy 
is mainly a consequence of the adiabatic app roximation. 

From the SOHO/GOLF seismic data (jBaudin et al.l 

120051) ■ we derive the maximum of the solar mode (intrinsic) 
surface velocity, that is 32.6 ± 2.6 cm/s. Then, using ^r, we 
infer the maximum of mode velocity at the photosphere , 
that is 18.5 ± 1.5 cm/s. According to lMichel etHI (|2009t) . 
the maximum of the solar mode (bolometric) amplitude in 
intensity is equal to 2.53 ±0.11 ppm. Then, by applying 
Eq. ([5]) in the case of the Sun, we derive the scaling factor 
f3 = 0.103 ± 10%. We have checked that this calibration 
depends very little on the choice of the chemical mixture 
(see also Sect. 14. 3p . We then adopt this value for the case 
of HD 49933. 



3.2. The mode intensity fluctuations measured by CoRoT 



The seismic analysis bv lBenomar et al.l (|2009[ ) provides the 
apparent amplitude Ai of the £ =0, 1 and 2 modes and 
the associated error bars. However, the CoRoT measure- 
ments Ai correspond to relative intensity fluctuations in 
the CoRoT passband. Furthermore, the observed (appar- 
ent) mode amplitudes depend on the degree £. Therefore, 
to transform them into bolometric and intrinsic intensity 
fluctuations normalised to the radial modes, we divide them 
by the CoRoT res ponse function, -R^, derived here for £ =0, 
1 and 2, following iMichel et all (|2009[ ). The adopted values 
for Ri are: Rq = 0.90, Ri = 1.10, and i?2 = 0.66. We finally 
derive the bolometric intensity fluctuations normalised to 
the radial modes according to: 



(SL/L) 



CoRoT 




We shall stress 
amplitudes derived by 



(7) 

tha t the differenc e s bet ween the 
Benomar et al.l (|2009D and by 



lAppourchaux et ahl (|2008D are smaller than the l-tr error 
bars. Furtherm ore, these amplitude s are in agreement with 
those found bv lMichel et al.l (|2008l ). using a different tech- 
nique. 

3.3. Comparison witli tlie CoRoT measurements 

We compute the mode amplitudes in terms of bolomet- 
ric intensity fluctuations, 5L/ L, according to Eqs. (O and 
([6]) (sec Sect. 13. 1|) . As for d^, the uncertainty associated 
with the measured mode linewidths, F, put uncertain- 
ties on the theoretical values of 5L/L. Furthermore, the 
uncertainty associated with the calibrated factor (3 (see 
Sect. 13. ip also puts an additional uncertainty on 5L/L. 
From here on, A{6L/L) will refer to the 1-a uncertainties 



((5i/L)y^(iAF/F) +(A/3//3)', where AF (reps. A/3) is 

the 1-a uncertainty associated with F (resp. /3). 

FigureHjcompares, as a function of the mode frequency, 
SL/L to the CoRoT measurements: {SL/L)coRoT- The dif- 
ference between our calculations and the observations is 
shown in the bottom panel. As for the velocity, this differ- 
ence must be compared with Ci, the 1-a interval resulting 
from the association of the 1-a error bars A{dL/L) and 
the 1-a error, A{6L/L)coRoT, associated with the CoRoT 
measurements. Accordingly, we have aL = y/a^ + W where 
a = A{5L/L) and 6 = A((5L/i)coRoT- 

As seen in Fig. [51 below v < 1.9 mHz, values of 
5L/L are within approximately 1-ctl in agreement with 
((5L/L)coRoT- However, above v ~ 1.9 mHz, the differences 
between 5L/L and ((5£/L)coRoT exceed 2-aL- 

Assuming a solar abundance ([Fe/H]= 0) results in a 
clear overestimation of A (^L/L)coRoT- Furthermo re, calcu- 
lations which assume the lGrevesse fc Noels! (jl993l) chemical 
mixture result in mode amplitudes larger by ~ 15%. 

Both in terms of intensity and velocity, differences 
between the calculated mode amplitudes and those de- 
rived from the observations (CoRoT and HARPS) are ap- 
proximately within the 1-a domain below i/ ~ 1.9 mHz. 
This then validates the intensity-velocity relation given by 
Eq. (O at the level of the current seismic precision . 

The maximum (SL/L) peaks at z/max — 1-9 mHz and 
the maximum of Vs at z^max — 1.8 mHz. By comparison, 
((5L/L)coRoT peaks at i^max — 1-8 mHz and uharps peaks 



at i/„ 



1.7 mHz. The difference in f,. 



between the 



observations (CoRoT and HARPS) and the model can be 
partially a consequence of the clear tendency at high fre- 
quency toward over-estimated amplitudes compared to the 
observations. 



4. Discussion 

4.1. Uncertainties in the knowledge of the fundamental 
parameters of HD 49933 

Uncertainties in the knowledge of Teff and log g place un- 
certainties on the theoretical values of T and hence on the 
mode amplitudes {vs and 5L/L). However, estimating these 
uncertainties would require the consideration of 3D mod- 
els with a Tcff and a log 5 that depart more than 1-a from 
the values adopted in our modehng, i.e. T^e = 6750 K and 
log g = 4.25. This is beyond the scope of our efforts since 
such 3D models are not yet available. 

4.2. Influence of the mode mass 



As discussed in details in I Samadi et al.l (j2008af ) , the com- 
puted mode surface velocities Vg significantly depend on the 
choice of the height h in the atmos phere where the mod e 
masses arc evaluated. According to ISamadi et al.l (l2008a[ ) , 
seismic measurements performed with the HARPS spectro- 
graph reflect conditions slightly below the formation depth 
of the K line. Accordingly, we have evaluated by default 
the mode masses at the optical depth where the K line is 
expected to be formed (i.e. T500 nm ~ 0.013), which corre- 
sponds, for our 3D models, to a height of about 350 km 
above the photosphere. We can evaluate how sensitive we 
are to the choice of h. Indeed, evaluating the mode mass 
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Fig. 2. Top: Mode bolomctric amplitude in intensity as a 
function of the mode frequency {v). The filled circles con- 
nected by the thick solid line correspond to the mode am- 
plitudes in intensity, 5L/L, derived for HD 49933 accord- 
ing to Eq. (O and Eq. ([T|) where the mode surface veloc- 
ity V is evaluated at the photosphere. The thick dashed 
solid line corresponds to the mode amplitude in intensity 
associated with the model with [Fe/H]= 0. The red tri- 
angles and associated error bars correspond to the mode 
amplitudes in inte nsity, {5L / L)coKoT, obtained by from 
the CoRoT data H enomar et all I2OO90 . These measure- 
ments have been translated into bolometric amplitudes fol- 
lowing |Micheretin] (|2009t ). Bottom: Same as top for the 
difference between 5L/L and ((5i/L)coRoT- The l-u error 
bars correspond here to y/a?~vT? where a = A{SL/L) and 
b = A((5L/L)coRoT (see text). 



at the photosphere results in values of Vs which are about 
15 % lower and hence would reduce the discrepancy with 
the HARPS observations. On the other hand, evaluating 
the mode mass one pressure scale height (~ 300 km at the 
photosphere) above h = 350 km results in an increase of Vs 
of about 10%. A more rigorous approach to derive the dif- 
ferent heights in the atmosphere where the measurements 
are sensiti ve would require a de dicated modeling (see a dis- 
cussion m ISamadi et aLll2008a[ ). 



4.3. The intensity-velocity relation 



Sensitivity to the location: 

The derivation of Eq. ^ (or equivalently Eq. ([5])) is based 
on the assumption that STcs cx ST\rp^rj^ (see Sect. 13. ip . 
This is quite a arbitrary simplification. In order to check 
how sensitive our results are to this assumption, we have 
computed Eq. (O and Eq. ^ for two different positions 
in the atmosphere. The first position, h = hi, is chosen 
one pressure scale height (~ 300 km) above the photo- 
sphere, which corresponds to an optical depth of t ~ 0.02. 
The second position, /i = /12, is chosen one pressure scale 
height beneath the photosphere, that is around r ~ 200. 
For both positions, the mode amplitudes with frequencies 
below ~ 1.9 mHz are almost unchanged. Concerning the 
amplitudes of modes with frequencies above ~ 1.9 mHz, 
they are increased by up to ~ 20 % when h = hi and are in 
turn almost unchanged when h = h2. Since the fluctuations 
of L induced by the oscillations are mostly due to temper- 
ature changes that occur around an optical depth of the 
order of the unit, we can conclude that our calculations are 
almost insensitive to the choice of the layer in the visible 
atmosphere where 6T is evaluated. 



Non-adiabatic effects: 

The modes are measured at the surface of the star where 
non-adiabatic interactions between the modes and convec- 
tion as well as radiative losses of the modes are important. 
Assuming Eq. ^ is then a crude approximation. In fact, 
it is clearly non-valid in the case of the Sun since it re- 
sults in a severe over-estimation of the solar mode ampli- 
tudes in intensity (see Sect. 13. 1| ). Avoiding this approxima- 
tion requires non-adiabatic eigenfunctions computed with 
a tin ie-dependent convection model. However, such models 
(e.g. iGrigahcene et all l2005l : lBalmfort"hl I1992D are subject 
to large uncertainties, and there is currently no consensus 
about the non-adiabatic mechanisms that play a si gnifi- 
cant role (see e.g. the recent review bv lHoudek|[2008[) . For 
instance, parameters are usually introduced in the theories 
so that they cannot be used in a predictive way. 

In the present study, we adopt by default the adiabatic 
approximation and introduce in Eq. ([5]) the parameter (3 
calibrated with helioseismic data. We show here that de- 
spite the deficiency of the quasi-adiabatic approximation, 
it nevertheless provides the correct scaling, at least at low 
frequency and at the level of the present seismic precisions. 

As an alternative approach, comparing the spectrum ob- 
tained from the 3D models in intensity with that obtained 
in velocity can provide valuable information concerning the 
intensity-velocity relation, in particular concerning the de- 
parture from the adiabatic approximation and the sensi- 
tivity to the surface metal abundance. We have started to 
carry out such a study. For the velocity, the (few) acous- 
tic modes trapped in the simulated boxes can be extracted 
and their properties measured. But this was impossible to 
do for the intensity with the simulations at our disposal 
because the computed spectrum for the intensity is domi- 
nated by the granulation background. As a consequence it 
is not possible to extract the mode amplitudes in intensity 
with sufficient accuracy . A comparison between the spec- 
tra obtained from the 3D models requires a much longer 
time series (work in progress). 
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Sensitivity to the metal abundance: 

We have shown in this study how the mode amplitudes 
in the velocity are sensitive to the surface metal abun- 
dance. An open question is how sensitive is the intensity- 
velocity relation in general to the metal abundance? A the- 
oretical answer to this question would require a realistic 
and validated non-adiabatic treatment. The pure numeri- 
cal approach mentionned above can also in principle pro- 
vide some answers to this question. However, as discussed 
above, this approach is not applicable with the time se- 
ries at our disposal. Concerning the quasi-adiabatic rela- 
tion of Eq. ([5]): a change of the metal abundance has a 
direct effect on Fa and an indirect effect on the properties 
of the (radial) eigen-displacement ^r- However, the com- 
parison between the metal-poor 3D model (SI) and the 3D 
model with the solar abundances (SO) shows that - at a 
fixed frequency Wosc ~ the ratio {SL/L)/v, which is equal 
to 4/3 (Fa — 1) (d^r/dr)/(^i. Wosc), is almost unchanged be- 
tween SO and SI (the differences are less than ^ 1%). In 
conclusion, the quasi-adiabatic relation of Eq. ^ depends 
very weakly on the surface metal abundance. Accordingly, 
the choice of the solar chemical mixture has a negligible 
impact on the value of the calibration factor /?. 



4.4. The solar case 

As seen in this study, the surface metal abundance has a 
pronounced effect on the mode excitation rates. One may 
then wonder about the previous validation of the theoret- 
ical model of stochast i c excitation in the ca se of the Sun 
(jBelkacem et all 120061: fSamadi et aLll2008bf) . Indeed, this 
validation was carried out with the use of a solar 3D model 
based on an "old" solar chemical mixt ure (namely those 
proposed by lAnders &: Grevess^ Il989l) w hile the "new" 
chemical mixture bv lAsplund et al.l (|2005f ) is characterized 
by a significantly lower metal abundance. 

In order to adress this issue, we have first considered two 
global ID solar models. One m odel has an " old" solar abun- 
dance (jGrevesse fc Noelj|l993l model Mqm hereafter) while 
the s econd one has the "new" abundances (jAsplund et al.l 
l2005l model Mnew hereafter). At the surface where the ex- 
citation occurs, the density of the solar model Mnow is only 
~ 5 % lower compared to the model Moid- According to 
the arguments developed in Paper I, this difference in the 
density must imply a difference in the convectivc velocities 
(u) of the order of-- (poid/pnow)^/^, where poid (resp. pncw) 
is the surface density associated with Mqm (resp. Mnew)- 
Accordingly, u is expected to be ~ 1.7 % higher for Mnew 
compared to Moid- 



The next question is what is the change in the solar 
mode excitation rates induced by the above difference in u7 
We have computed the solar mode excitation rates exactly 
in the same manner as for HD 49333 by using a solar 3D 
simulations based on the " old" abundances. We obtained a 
rather good agre ement with th e different helioseismic data 
(see the result in lSamadi et al. 20 08b. ) ■ To derive the solar 
mode excitations expected with the "new" solar abudance, 
we have proceeded in a similar way as the one done in Paper 
I: we have increased the convective velocity u derived from 
the solar 3D model by 2 % while keeping the kinetic flux 
constant (see details in Paper I). This increase of ^ 2 % of 
u results in an increase of ~ 10 % of the mode excitation 
rates. This increase is significantly lower than the current 



uncertainties associated with the different helioseismic data 
(jBaudin et al.ll2005l : ISamadi et ani2008bD . 

4.5. Discrepancy at high frequency 

The discrepancy betwen theoretical calculations and ob- 
servations is particularly pronounced at high frequency. 
This discrepancy may be attributed to a canceling be- 
tween the entropy and the Reynolds stress contributions 
(see Sect. I4.5.ip or the "scale length separation" assump- 
tion (see Sect. 14. 5. 2p . 

4.5.1. Canceling between the entropy and the Reynolds stress 
contributions 

The relative contribution of the entropy fluctuations to the 
excitation is found to be about 30 % of the total excitation. 
This is two times larger than in the case of the Sun (~15 %). 
This can be explained by the fact HD 49933 i s significantly 
hotter than the Sun and, as pointed-out by I Samadi et al.l 
(|2007[ ). the larger (L/M) oc T^g/ff, the more important the 
relative contribution of the entropy. Although more impor- 
tant than in the Sun, the contribution of the entropy fluctu- 
ations remains relatively smaller than the uncertainties as- 
sociated with the current seismic data. This is illustrated in 
Fig. [3J the difference between theoretical mode amplitudes 
which take into account only the Reynolds stress contribu- 
tion (C^, see Eq. (3) of Paper I) and those that include 
both contributions (entropy and Reynolds stress) is lower 
than ffy. In terms of amplitudes, the entropy fluctuations 
contribute only ~ 15 % of the global amplitude. This is 
significantly smaller than the uncertainties associated with 
the current seismic measurements. Seismic data of a better 
quality are then needed to constrain the entropy contribu- 
tion and its possible canceling with the Reynolds stress. 

Numerical simulations show some cancellation between 
the en tropy source term and the one due to the Reynolds 
stress (jStein et al.ll2004) . However, in the present theoreti- 
cal model of stochastic excitation, the cross terms between 
the entropy fluctuations and the Reynolds stresses vanish 
fsee ISamadi fc Goupill [200ll) . This is a consequence of the 
diffe rent assumptions conce rning the entropy fluctuations 
fs ee ISamadi fc Goupilll200ll see also the recent discussion 
in lSamadi et al.l ( 2008H )). Accordinglv. the entropy source 
term is included as a source indepen dent frorn the R eynolds 
stress contribution. As suggested bv lHoudeS (|2006( ). a par- 
tial canceling between the entropy fluctuations and the 
Reynolds stress can decrease the mode amplitudes of F- 
type stars and reduce the discrepancy between the theoret- 
ical calculations and the observations. 

There is currently no theoretical description of these in- 
terferences. In order to have an upper limit of the interfer- 
ences, we assume that both contributions locally and fully 
intcrfer. This assumption leads to the computation of the 
excitation rates per unit mass as: 



dm 



dy 

dm 



dy 

dm 



where {d'y/dm)p,g and (dT/dm)^ arc the contributions 
per unit mass of the Reynolds stress and entropy respec- 
tively. The result is presented in Fig. |3| in terms of velocity 
(top panncl) and in terms of intensity (bottom pannel). 
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The mode amplitudes are decreased by up to '-^ 55%. In 
that case, ((5L/L)coRoT is systematically under-estimated. 
Obviously, a partial canceling between the entropy contri- 
bution and the Reynolds stress would result in a smaller 
decrease. 

We have assumed here that the cancellation between 
the two terms is independent of the mode frequ ency (see 
Eq. However, according to lStein et~aLl (|2004[ ). the level 
of the cancellation depends on the frequency (see their 
Fig. 8). In particular, for F-type stars, the cancellation is 
expected to be more important around and above the peak 
frequency. 

As a conclusion, the existence of a partial canceling 
between the entropy fluctuations and the Reynolds stress 
can decrease the mode amplitude and could improve the 
agreement with the seismic observations at high frequency. 
However, there is currently no theoretical modeling of the 
interference between theses two terms. Further theoretical 
developements are required. 

4.5.2. The "scale length separation" assumption 

Th e "scale length separ ation" assumption (see the review 
by I Samadi et al.l l2008bf ) consists of the assumption that 
the eddies contributing effectively to the driving have a 
characteristic length scale smaller than the mode wave- 
length. This assumption is justified for a low Mach num- 
ber (Aft). However, this approximation is less valid in the 
super-adiabatic region where Alt reaches a maximum (for 
the Sun Mt is up to 0.3) and accordingly affects the high- 
frequency modes more. This approximation is then ex- 
pected to be even more questionable for stars hotter than 
the Sun, since Mt increases with Tcff. This spatial separa- 
tion can be avoided, however if the kinetic energy spectrum 
associated with the turbulent elements {E{k)) is properly 
coupled with the spatial dependence of the modes (work 
in progress). In that case, we expect a more rapid decrease 
of the driving efficiency with increasing frequency than in 
the present formalism where the spatial dependence of the 
modes is totally decoupled from E{k) (i.e. "scale length 
separation"). 

5. Conclusion 

From the mode linewidths measured by CoRoT and theo- 
retical mode excitation rates derived for HD 49933, we have 
derived the expected mode surface velocities Vg which we 
have compared with wharps, the mode velocity spectrum 
derived from the seis mic observations ob tained with the 
HARPS spectrograph (|Mosser et al.ll2005[ ). Except at high 
frequency > 1.9 mHz), the agreement between computed 
Vs and Wharfs is within the 1-ct domain associated with 
the seismic data from the HARPS spectrograph. However, 
there is a clear tendency to overestimate wharfs above 
v ^ 1.9 mHz,. 

Using a calibrated quasi-adiabatic approximation to re- 
late the mode velocity to the mode amplitude in inten- 
sity (Eq. O, we have derived for the case of HD 49933 the 
expected mode amplitudes in intensity. Computed mode 
intensity fluctuations, SL/L, are within 1-ct in agreement 
with the seismic cons traints derived from the CoRoT data 
(jBenomar et al.|[2009f ). However, as for the velocity, there is 
a clear tendency at high frequency {v ^ 1.9 mHz) towards 
over-estimated SL/L compared to the CoRoT observations. 





Fig. 3. Top: Same as Fig. [TJ The thin dashed line corre- 
sponds to a calculation that takes only the contribution of 
the Reynolds stress into account. The dot-dashed line corre- 
sponds to a calculation in which we have assumed that the 
contribution of the Reynolds stress interferes totally with 
that of the entropy fluctuations (see text). The thick solid 
line has the same meaning as in Fig[T] Bottom: Same as 
top for 6L/L. The triangles and associated error bars have 
the same meaning as in Fig. [2] 



Calculations that assume a solar surface metal abun- 
dance result, both in velocity and in intensity, in ampli- 
tudes larger by ~ 35 % around the peak frequency (^'max — 
1.8 mHz) and by up to a factor of two at lower frequency. 
It follows that, ignoring the current surface metal abun- 
dance of the star results in a more severe over-estimation 
of the computed amplitudes compared with observations. 
This illustrates the importance of taking the surface metal 
abundance of the solar-like pulsators into account when 
mo deling the mode driving. In addition, we point out that 
the .Grevesse fc Noelj (|l993f ) solar chemical mixture results 
in mode amplitudes larger by about 15 % with respect to 
calculations that ass ume the "new" solar abundance by 
lAsplund et al.l (j2005l ). However, this increase remains sig- 
nificantly smaller than the uncertainties associated with 
current seismic measurements. 

Since both mode amplitudes in terms of surface velocity 
and intensity are available for this star, it was possible to 
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test the vahdity of the cahbrated quasi-adiabatic relation 
(Eq. (O). Our comparison shows that this relation provides 
the correct scaling, at least at the level of the present seis- 
mic precisions, . 

Both in terms of surface velocity and of intensity, the 
differences between predicted and observed mode ampli- 
tudes are within the 1-a uncertainty domain, except at 
high frequency. This result then validates for low frequency 
modes the basic underlying physical assumptions included 
in the theoretical model of stochastic excitation for a star 
significantly different in effective temperature, surface grav- 
ity, turbulent Mach number (Mt) and metallicity compared 
to the Sun or a Cen A. 

As discussed in Sect. SI the clear discrepancy between 
predicted and observed mode amplitudes seen at high fre- 
quency may have two possible origins: First, a canceling 
between the entropy contribution and the Reynolds stress 
is expected to occur and to be important around and above 
the frequency of the maximum of the mode excitation rates 
(see Sect. I4.5TT|) . Se cond, the assumption called the "scale 
length separation" (ISamadi et all l2008bD may also result 
in an over-estimation of the mode amplitudes at high fre- 
quency (see Sect. 14.5."^ . These issues wiU be investigated 
in a forthcoming paper. 
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